      SUBROUTINE SETUP
      RETURN
      END
C   HS9
      BLOCK DATA
      INCLUDE 'BLOCO.INC'
      INTEGER I,J
      DATA NAME/'HS9 '/
      DATA X/  NX*0.0D0/
      DATA N/2/ , NH/1/ , NG/ 0/
      DATA DEL0/1.00D0/ ,TAU0/1.D0/ ,TAU/.1D0/
      DATA ((GUNIT(I,J),I=1,3),J=0,2)/-1,0,0,-1,0,0,-1,0,0/
      END
      SUBROUTINE EF(X,FX)
      INCLUDE 'FUCO.INC'
      DOUBLE PRECISION X(NX),PI,FX
      SAVE PI
      DATA PI/3.141592653587D0/
      ICF=ICF+1
      FX=DSIN(PI*X(1)/12.D0)*DCOS(PI*X(2)/16.D0)
      RETURN
      END
      SUBROUTINE EGRADF(X,GRADF)
      INCLUDE 'FUCO.INC'
      DOUBLE PRECISION X(NX),GRADF(NX)
      DOUBLE PRECISION  PI
      SAVE PI
      DATA PI/3.141592653587D0/
      ICGF=ICGF+1
      GRADF(1)=DCOS(PI*X(2)/16.D0)*PI/12.D0*DCOS(PI*X(1)/12.D0)
      GRADF(2)=-DSIN(PI*X(1)/12.D0)*PI/16.D0*DSIN(PI*X(2)/16.D0)
      RETURN
      END
      SUBROUTINE EH(I,X,HXI)
      INCLUDE 'FUCO.INC'
      INTEGER I
      DOUBLE PRECISION X(NX),HXI
      CRES(I)=CRES(I)+1
      HXI= 4.D0*X(1)-3.D0*X(2)
      RETURN
      END
      SUBROUTINE EGRADH(I,X,GRADHI)
      INCLUDE 'FUCO.INC'
      INTEGER I
      DOUBLE PRECISION X(NX),GRADHI(NX)
      CGRES(I)=CGRES(I)+1
      GRADHI(1)=4.D0
      GRADHI(2)=-3.D0
      RETURN
      END
      SUBROUTINE EG(I,X,GXI)
      INCLUDE 'FUCO.INC'
      INTEGER I
      DOUBLE PRECISION X(NX),GXI
      RETURN
      END
      SUBROUTINE EGRADG(I,X,GRADGI)
      INCLUDE 'FUCO.INC'
      INTEGER I
      DOUBLE PRECISION X(NX) ,GRADGI(NX)
      RETURN
      END


